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Measurements of the mass and width of the rjc using ^/'(3686) — ?■ 7?7c 
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The mass and width of the lowest-lying 5- wave spin singlet charmonium state, the r^c, are mea- 
sured using a data sample of 1.06 x 10* ^(3686) decays collected with the BESIII detector at the 
BEPCII storage ring. We use a model that incorporates full interference between the signal reac- 
tion, ^(3686) — ^ 7»7c, and a non-resonant radiative background to describe the line shape of the 77c 
successfully. We measure the 77c mass to be 2984.3 ± 0.6 ± 0.6 MeV/c^ and the total width to be 
32.0 ± 1.2 ± 1.0 MeV, where the first errors are statistical and the second are systematic. 

PACS numbers: 13.25.Gv, 13.20.Gd, 14.40.Pq 



In recent years, many new charmonium or 
charmonium-like states have been discovered. These 
states have led to a revived interest in improving the 
quark- model picture of hadrons [ij. Even with these 
new discoveries, the mass and width of the lowest-lying 
charmonium state, the 77c, continue to have large uncer- 
tainties when compared to those of other charmonium 
states Early measurements of the properties of 

the Tjc using J/ij; radiative transitions [3, |j| found 
a mass and width near 2978 MeV/c^ and 10 MeV, 
respectively. However, recent experiments, including 
photon-photon fusion and B decays, have reported a 
significantly higher mass and a much larger width 0- 3 . 
The most recent study by the CLEO-c experiment 9[, 
using both ■0(3686) — 777c and J/t/i — 777c, pointed 
out a distortion of the r?c line shape in -0(3686) decays. 
CLEO-c attributed the r]c line-shape distortion to the 
energy dependence of the Ml transition matrix element. 

In this Letter, we report measurements of the rjc mass 



and width using the radiative transition ■0(3686) — >■ 777c. 
We successfully describe the measured rjc line shapes 
using a combination of the energy dependence of the 
hindered-Ml transition matrix element and a full inter- 
ference with non-resonant ^0(3686) radiative decays. The 
analysis is based on a 0(3686) data sample of 1.06 x 10* 
events [l3| collected with the BESIII detector operating 
at the BEPCII e+e" collider. A 42 pb~^ continuum data 
sample, taken at a center-of-mass energy of 3.65 GeV, is 
used to measure non-^0(3686) backgrounds. 

The mass and width are determined from fits to 
the invariant mass spectra of exclusive 77c decay modes. 
Six modes are used to reconstruct the 77c: KsK'^tt~ , 

and 3(7r^"7r^), where the Ks is reconstructed in tt+tt^, 
and the 77 and tt*^ in 77 decays. The inclusion of charge 
conjugate modes is implied. 

The BESIII detector is described in detail in Ref. [ill ■ 

The detector has a geometrical acceptance of 93% of 
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4tt. a small cell helium-based main drift chamber 
(MDC) provides momentum measurements of charged 
particles; in a 1 T magnetic field the resolution is 0.5% at 
1 GeV/c. It also supplies an energy loss {—dE/dx) mea- 
surement with a resolution better than 6% for electrons 
from Bhabha scattering. The electromagnetic calorime- 
ter (EMC) measures photon energies with a resolution 
of 2.5% (5%) at 1 GeV in the barrel (endcaps). The 
time-of-flight system (TOF) is composed of plastic scin- 
tillators with a time resolution of 80 ps (110 ps) in the 
barrel (endcap) and is mainly useful for particle identi- 
fication. The muon system provides 2 cm position reso- 
lution and measures muon tracks with momenta greater 
than 0.5 GeV/c. 

We use inclusive Monte Carlo (MC) simulated events 
as an aid in our background studies. The V'(3686) reso- 
nance is produced by the event generator KKMC [l^ . 
while the decays are generated by EvtCen |13| with 
known branching fractions or by Lundcharm for 
unmeasured decays. The signal is generated with an an- 
gular distribution of 1 -I- cos^ 9^ for ■(/;(3686) ^rjc, and 
phase space for multi-body rjc decays, where 9^ is the 
angle between the photon and the positron beam direc- 
tion in the center-of-mass swtem. Simulated events are 
processed using GEANT4 15|, where measured detector 
resolutions are incorporated. 

We require that each charged track (except those from 
Ks decays) is consistent with originating from within 
1 cm in the radial direction and 10 cm along the beam 
direction of the run- by-run-determined interaction point. 
The tracks must be within the MDC fiducial volume, 
I cos6'| < 0.93. Information from the TOF and —dE/dx 
is combined to form a likelihood C^^ (or Ck) for a pion (or 
kaon) hypothesis. To identify a track as a pion (kaon), 
the likelihood C^^ {^k) is required to be greater than 
0.1% and greater than C)c [Ct^). 

Photons are reconstructed from isolated showers in the 
EMC that are at least 20 degrees away from charged 
tracks. The energy deposited in the nearby TOF scin- 
tillator is included to improve the reconstruction effi- 
ciency and the energy resolution. Photon energies are 
required to be greater than 25 MeV in the fiducial EMC 
barrel region (|cos0| < 0.8) and 50 MeV in the end- 
cap (0.86 < |cos0| < 0.92). The showers close to the 
boundary are poorly reconstructed and excluded from 
the analysis. Moreover, the EMC timing, with respect 
to the collision, of the photon candidate must be in co- 
incidence with collision events, i.e. < t < 700 ns, to 
suppress electronic noise and energy deposits unrelated 
to the event. 

The Ks — >■ 7r+7r~ candidates are reconstructed from 
pairs of oppositely charged tracks. The secondary ver- 
tex constrained tracks must have an invariant mass 
±10 MeV/c^ of the nominal Ks mass, and a decay length 
more than twice the vertex resolution. The track in- 



formation at the secondary vertex is an input to the 
kinematic fit. Candidate tt° and rj mesons are recon- 
structed from pairs of photons with an invariant mass in 
the range 0.118 GeV/c^ < M{-/-f) < 0.150 GeV/c^ for 
TfO and 0.50 GeV/c^ < Mijj) < 0.58 GeV/c^ for r?. The 
remaining photons are considered as candidates of the 
transition photon. 

Events with either extra charged tracks or non-zero net 
charge are rejected. The rjc candidates are reconstructed 
from KsK^TT^ , K^K^tt^, rjir^Tr^, KsK^tt^tt^ tt^ , 
K^K^TT^TT~TT^ and 3(7r^7r~). We select events in the 
region 2.7 GeV/c^ < M{r]c) < 3.2 GeV/c^. A four- 
constraint (4C) kinematic fit of all selected charged parti- 
cles and the transition photon with respect to the initial 
'0(3686) four-momentum is performed to reduce back- 
ground and improve the mass resolution. When addi- 
tional photons are found in an event, we loop over all 
possible combinations and keep the one with the best 
xlc from the kinematic fit. The xic required to be 
less than 60, a value is determined by optimizing the fig- 
ure of merit for most of the channels, S/\/S + B, where 
S {B) is the number of signal (background) events in the 
signal region (2.9 GeV/c^ < M{r]^) < 3.05 GeV/c^). In 
addition, to remove '0(3686) — >■ tt^tt~ J/tp events, we re- 
quire there be no tt+tt^ pair with a recoil mass in the J/tp 
signal region. To suppress background from tt*^ —^77, we 
demand that the transition photon should not form a 7r° 
with any other photon in the event. 

The main source of background is from '0(3686) — 
Tr'^Xi decays, where a photon from the tt*^ — > 77 de- 
cay is missing, and Xi represents the r]c final states un- 
der study. These decays could proceed via various in- 
termediate states, and most of the branching fractions 
are unknown. To estimate their contribution, we recon- 
struct '0(3686) — >■ TT^Xi decays from data. The selec- 
tion criteria are similar to those applied to the 777c can- 
didates except an additional photon is required. The 
'0(3686) — TT^Xi signal yields are extracted from fits 
to the M(77) invariant mass distributions for different 
M{Xi) mass bins. The relative efiiciencies of the jrjc and 
TT^Xi selection criteria are estimated in each M (Xi) mass 
bin using phase space distributed 0^(3686) — > ir^Xi MC 
events. Combining this relative efficiency with the num- 
ber of '0(3686) — ?> Tr'^Xi signal events in every M{Xi) bin, 
we estimate the n^Xi events that pass the 777c selection. 
We also examine the efficiencies ofn'^Xi events generated 
with different dynamics, and the change is negligible. 

Other potential 0'(3686) decay backgrounds are in- 
vestigated using 1.06 X 10* inclusive MC events where 
TT^Xi events have been excluded. We find no other dom- 
inant background processes, but do find dozens of de- 
cay modes that each makes small additional contribu- 
tions to the background. These decays typically have 
additional or fewer photons in their final states. The 
sum of these background events is used to estimate the 
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contribution from other ■!/)(3686) decays. Backgrounds 
from the e^e~ — > qq continuum process are studied us- 
ing a data sample taken at -/s = 3.65 GeV. Contin- 
uum backgrounds are found to be small and uniformly 
distributed in M{Xi). There is also an irreducible non- 
resonant background, ?/'(3686) — > 7^^, that has the same 
final state as signal events. A non-resonant component 
is included in the fit to the rjc invariant mass. 

Figure [1] shows the rjc invariant mass distributions for 
selected rjc candidates, together with the estimated tt^Xi 
backgrounds, the continuum backgrounds normalized by 
luminosity, and other ?/'(3686) decay backgrounds esti- 
mated from the inclusive MC sample. A clear rjc signal 
is evident in every decay mode. We note that all of the 
•qc signals have an obviously asymmetric shape: there is 
a long tail on the low-mass side; while on the high-mass 
side, the signal drops rapidly and the data dips below 
the expected level of the smooth background. This be- 
havior of the signal suggests possible interference with 
the non-resonant ^Xi amplitude. In this analysis, we as- 
sume 100% of the non-resonant amplitude interferes with 
the Tjc- 

The solid curves in Fig. [1] show the results of an un- 
binned simultaneous maximum likelihood fit in the range 
from 2.7 to 3.2 GeV/c^ with three components: signal, 
non-resonant background, and a combined background 
consisting oiir'^Xi decays, continuum, and other ?/'(3686) 
decays. The signal is described by a Breit-Wigner func- 
tion {BW) convolved with a resolution function. The 
non-resonant amplitude is real, and is described by an 
expansion to second order in Chebychev polynomials de- 
fined and normalized over the fitting range. The com- 
bined background is fixed at its expected intensity, as 
described earlier. The fitting probability density func- 
tion (PDF) as a function of mass (to) reads: 



F{m) 



e(m) 



TO 



6 (to) 



where S{m), J\f{m) and B{m) are the signal, the non- 
resonant jXi component, and the combined background, 
respectively; E-y is photon energy; a is the experimen- 
tal resolution and e(TO) is the mass-dependent efficiency. 
The EJ^ multiplying |5(to)P reflects the expected energy 
dependence of the hindered-Ml transition 16], which 
partially contributes to the 77c low mass tail as well as 
the interference effect. The interference phase, (j), and 
the strength of the non-resonant component, a, are al- 
lowed to vary in the fit. 

The mass-dependent efficiencies are determined from 
phase space distributed MC simulations of the 77c decays. 
Efficiencies obtained from MC samples that include in- 
termediate states change the resulting mass and width 
by negligible amounts. MC studies indicate that the res- 
olution is almost constant over the fitting range. Thus, 
a mass-independent resolution is used in the fit. The 



detector resolution is primarily determined by MC sim- 
ulation for each r]c decay mode. The consistency be- 
tween the data and MC simulation is checked by the de- 
cay ■0(3686) — >■ 77J/'0, where the J/tp decays into the 
same final states as the rjc- We use a smearing Gaussian 
function to describe possible discrepancies between data 
and MC simulations. By fitting the MC-determined J/t/j 
shape convolved by a smearing Gaussian function to the 
data, we determine the parameters of the Gaussian func- 
tion. Due to the different kinematics, the parameters are 
slightly different for each mode. 

In the simultaneous fit, the rjc mass and width are con- 
strained to be the same for all the decay modes but still 
free parameters; the two Chebyshev polynomial coeffi- 
cients and the factor a are also allowed to float. Two 
solutions for the relative phase are found for each decay 
mode, one corresponds to constructive and the other de- 
structive interference between the two amplitudes at the 
77c peak. Regardless of which solution we take, the mass, 
width of the 77^ and the overall fit quality are always un- 
changed [111. The mass is M = 2984.3 ± 0.6 MeV/c^, 
and width F = 32.0 ± 1.2 MeV. The goodness-of-fit 
X^/ndf — 283.4/274, which indicates a reasonable fit. 
The solutions for relative phase of each mode are listed 
in Table m 

TABLE I: Solutions of relative phase (in unit of radian) of 
each decay mode. 



mode 


constructive 


destructive 


KsK+n- 


2.94 ±0.27 


3.75 ±0.26 




2.63 ±0.21 


3.96 ±0.19 




2.41 ±0.13 


4.28 ±0.09 


K S K'^ Tt^ ■K~ 'R~ 


2.16 ±0.11 


4.46 ±0.07 


K+K-n+TT--!v° 


2.73 ±0.19 


4.00 ±0.16 


3{-K + TV-) 


2.28 ±0.10 


4.43 ± 0.06 



However, without the interference term, the fit would 
miss some data points, especially where the symmetric 
shape of a Breit-Wigner function is deformed, and the 
goodness-of-fit is /ndf = 426.6/280. The statistical 
significance of the interference, calculated based on the 
differences of likelihood and degrees of freedom between 
fits with and without interference, is of order 15(7. 

The systematic uncertainties of the 77c mass and width 
mainly come from the background estimation, the mass 
scale and resolution, the shape of the non-resonant com- 
ponent, the fitting range, and the efficiency. 

In the fit, the ir^Xi background is fixed at its expected 
intensity, so the statistical uncertainty of the observed 
Tr'^Xi events introduces a systematic error. To estimate 
this uncertainty, we vary the number of events in each 
bin by assuming Gaussian variations from the expected 
value. We repeat this procedure a thousand times, and 
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FIG. 1: The M{Xi) invariant mass distributions for the decays KsK^ ■k~ , K~ tv^ , rj-K^-K~ , KsK^-k'^tx~-k~ , K'^ K~ n^7v~ tv" 
and 3(71"*" 7r~), respectively, with the fit results (for the constructive solution) superimposed. Points are data and the various 
curves are the total fit results. Signals are shown as short-dashed lines, the non-resonant components as long-dashed lines, 
and the interference between them as dotted lines. Shaded histograms are (in red/yellow/green) for (continuum/vr'^Xi/other 
'i/>(3686) decays) backgrounds. The continuum backgrounds for KsK^-k^ and •qn'^-K' decays are negligible. 



take the standard deviation of the resulting mass, width, 
and phases as systematic errors. We also use different 
dynamics in generating the tt'^X^ events (with the same 
final state, but different intermediate states) for the ef- 
ficiency correction, and find the differences in resulting 
mass and width are small. We take 0.24 MeV/c^ in mass 
and 0.44 MeV in width as the systematic errors for the 
TT^Xi background estimation. 

We assign a 0.07 MeV/c^ (0.06 MeV) error in mass 
(width) for the non-resonant component shape that is 
obtained by changing the polynomial order. Also we in- 
clude an additional non- interfering component, which is 
represented by a 2nd-order polynomial with free strength 
and shape parameters. The changes in the resulting rjc 
mass and width are 0.10 MeV/c^ and 0.02 MeV, respec- 
tively; and the fraction of this component to total non- 
resonant rate varies from to 25% depending on decay 
mode. These variations are included in systematic errors. 

The systematic error from the uncertainty in the other 
V'(3686) decay backgrounds is estimated by floating the 
magnitude and changing the shape of this component 
to a 2nd-order polynomial with free parameters. The 
changes, 0.05 MeV/c^ in mass and 0.06 MeV in width, 
are taken as systematic errors. 

The consistency of the mass scale and resolution be- 
tween data and MC simulations is checked with the de- 
cay V'(3686) — > jjJ/ip, and possible discrepancies are 
described by a smearing Gaussian distribution, where a 



non-zero mean value indicates a mass offset, and a non- 
zero a represents difference between the data and MC 
mass resolutions {cr'^i^^^ — c^mc)^^^- ^ typical mass shift 
is about —1.0 MeV/c^ and resolution smear is ^ 3.0 
MeV. Another possible bias is the difference between 
input and the value after event-reconstruction and selec- 
tion. This is small for both the mass shift (< 0.3 MeV) 
and resolution smear. Both of these are added in the 
smearing Gaussian distribution. By varying the param- 
eters of the smearing Gaussian distribution from the ex- 
pected value, we estimate the uncertainties. From a large 
number of tests, the standard deviation of the resulting 
mass (width), 0.38 MeV/c^ (0.27 MeV), is taken as a 
systematic error in mass (width) for the mass scale un- 
certainty. A 0.35 MeV/c^ (0.60 MeV) systematic error 
in mass (width) is assigned due to the mass resolution 
uncertainty. 

The systematic error due to the fitting range is 
estimated by varying the lower-end between 2.6 and 
2.8 GeV/c^ and the higher-end between 3.1 and 
3.3 GeV/c^. The changes, 0.05 MeV/c^ in mass and 
0.07 MeV in width, are assigned as systematic errors. A 
mass-dependent efficiency is used in the fit. By remov- 
ing the efficiency correction from the fitting PDF, the 
changes, which are 0.05 MeV/c^ in mass and 0.06 MeV 
in width, are taken as systematic errors. The stability 
of the simultaneous fit program is checked by repeating 
the fit a thousand times with random initialization; the 
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standard deviation of mass and width, 0.14 MeV/c^ and 
0.66 MeV, respectively, are taken as systematic errors. 

We assume all these sources are independent and take 
their sum in quadrature as the total systematic error. We 
obtain the rjc mass and width to be 

M = 2984.3 ± 0.6 ± 0.6 MeV/c^ 

r = 32.0 ± 1.2 ± 1.0 MeV. 

Here (and elsewhere) the first errors are statistical and 
the second are systematic. 

The relative phases for constructive interference or de- 
structive interference from each mode are consistent with 
each other within Scr, which may suggest a common phase 
in all the modes under study. A fit with a common 
phase {i.e. the phases are constrained to be the same) 
describes the data well, with a y^/ndf = 303.2/279. 
Comparing to the fit with separately varying phases for 
each mode, we find the statistical significance for the 
case of five distinct phases to be 3.1cr. This fit yields 
M = 2983.9±0.6±0.6 MeV/c^, T = 31.3±1.2±0.9 MeV, 
and (/) = 2.40 ± 0.07 ± 0.47 rad (constructive) or = 
4.19 ± 0.03 ± 0.47 rad (destructive). The physics behind 
this possible common phase is yet to be understood. 

In summary, we measure the rjc mass and width via 
■0(3686) — ?> 777c by assuming all radiative non-resonant 
events interfere with the rjc. These results are so far the 
most precise single measurement of the mass and width 
of 77c ^l- For the first time, interference between the 
77c and the non-resonant amplitudes around the 77c mass 
is considered; given the assumptions of our fit, the sig- 
nificance of the interference is of order IScr. We note 
that this interference affects the rjc mass and width sig- 
nificantly, and may have impacted all of the previous 
measurements of the rjc mass and width that used ra- 
diative transitions. Our results are consistent with those 
from photon- photon fusion and B decays 041] ; this may 
partly clarify the discrepancy puzzle discussed above. 
The changes of the rjc mass and width may also have 
an impact on the expected 77^ mass and width, and will 
modify the parameters used in charmonium potential 
models, where the 77c mass is one of the input parame- 
ters. From this measurement, we determine the hyperfine 
mass splitting to be AM,,/(lS')cc = M{J/ip) - M{r]c) = 
which agrees well with recent lat- 
- 2^ as well as quark model predic- 
tions [21], and sheds light on spin-dependent interactions 
in quarkonium states. 
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